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Abstract — This paper presents a comprehensive non-linear 
dynamic model of a solid oxide fuel cell (SOFC) that can be 
used for transient behaviors studies. The model based on 
electrochemical and thermal equations, accounts for 
temperature dynamics and output voltage losses. The 
relaxation time is strongly related to the transient temperature 
distribution of the solid oxide fuel cell structure. Therefore, 
it is in the order of some minutes depending on the design 
parameters and the operating conditions. The model contains 
the hydrogen, oxygen and water block separately. Other blocks 
are concentration, activation and ohmic losses block. This 
analysis is based on an integrated dynamic model of the entire 
power plant using SIMULINK in Matlab. The analytical details 
of how active and reactive power output of a stand-alone solid 
oxide fuel cell power plant (FCPP) is controlled. 

Index Terms— Fuel cell, FCPP, SOFC. 

I. Introduction 

Many researchers reported on of Molten carbonate fuel 
cell. Looking back, in 2004, Francisco Jurada presented a 
model of 'Solid oxide fuel cell' without considering the 
temperature effect. The concept of molten carbonate fuel cell 
is exactly similar to solid oxide fuel cell. This comparison is 
given in [ 1 ] and [2] . The fuel input to the fuel cell is changing 
in steps in a paper in the literature [3], Another very useful 
model is available in the literature in which was suggested 
by Kourosh Sedghisigarchi, [3] and [4] .This paper considers 
all the subsystems of fuel cell, that includes hydrogen block, 
oxygen block, water block, activation and concentration block, 
and temperature block.. In temperature block, the relaxation 
time is closely related to the transient temperature distribution 
of the solid oxide fuel cell structure. The internal cell 
resistances are strongly temperature dependent. Thus, the 
relaxation time depends on the thermal properties, size and 
configuration of the cell, and operating conditions. In [5] and 
[3], SOFC model is given without considering thermal unit. 
However, in this paper, these structures are modified by the 
modeling thermal unit. A first comprehensive nonlinear 
dynamic model of solid oxide fuel cell that can be used for 
dynamic and transient stability studies is developed by K. 
Sedghisigarchi, Ali Feliachi in 2004 [6] . The model based on 
electrochemical and thermal equations, accounts for 
temperature dynamics and voltage losses. The output voltage 
response of a stand-alone fuel-cell plant to a step change, a 
fuel flow step change, and fast load variations are simulated 
to illustrate the dynamic behavior of SOFC for fast and slow 
perturbations. 

This paper presents a SOFC model and designs the control 
strategies for the AC voltage control and the active/reactive 
power control of the DC/AC inverter. Two separate 
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controllers are designed for these purposes. The fuzzy logic 
control scheme is employed for the design of the two 
controllers. 

II. Principles of fuel cell model 
The fuel cell (FC) is an electrochemical device that 
converts chemical energy of hydrogen gas ( H 2 ) and oxygen 

gas ( 2 ) into electrical energy. The solid oxide fuel cell 

consists of two porous ceramic electrodes separated by a 
dense ceramic electrolyte. The cell produces electricity by 
the electrochemical reaction of fuel (hydrogen and / or carbon 
monoxide) and the oxidant (oxygen) across the solid 
electrolyte. Oxygen fed to the air electrode (cathode) accepts 
electrons from the external circuit to form oxygen ions. The 
ions are conducted through the solid electrolyte to the fuel 
electrode (anode). At the fuel electrode, the ions combine 
with hydrogen and/ or carbon monoxide in the fuel to form 
water and /or carbon monoxide. This reaction releases 
electrons. Electrons flow from the fuel electrode (anode) 
through the external circuit back to the air electrode 
(cathode). The overall reaction is exothermic; the cell produces 
heat in addition to electricity. 

A typical fuel cell reaction is given below. The chemical 
reactions inside the cell that are directly involved in the 
production of electricity are given as 
At Anode: 

H 2 +0" -> H 2 + 2e~ 
CO + O '""->• CO 2 +2e" 
At Cathode: 



2 + 4e" -> 20 2 



Overall: 



H 2 +0 2 + CO^ H 2 + C0 2 




2H +2e +1/20,— »II,0 



Fig. 1 Basic operation of Fuel Cell 

The SOFC model consists of a) Electrochemical model- 
component material balance equations, b) Thermal model- 
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Energy balance equations, c) Nernst voltage equation d) 
Reformer model, e) Power conditioning unit model, f) Active 
and reactive power control. 

A. Component material balance equation 

The change in concentration of each species that appears 
in the chemical reaction can be written in terms of input and 

output flow rates and exit molarities (x i ) due to the following 

chemical reaction 

PV d 



N , 



iV ; ° +N/ 



(1) 



RT i dt 

where, V is compartment volume ( y^ ), N { , N { are molar 
flow rates (mol/s) of j* reactant at the cell input and output 

(exit), respectively. N- is the reaction rate (mol/s) of the 

reactant. is the cell temperature in degree Kelvin, is cell 
pressure (atom), is the gas constant [( 1 atom)/ (k mol degree 
Kelvin)]. 

Molar flow of any gas through the valve is proportional 
to its partial pressure inside the channel. The molar expression 
is given by 



N 



Hi 



K f 



N 



H 2 



K 



H 2 



where, N H ^ is the hydrogen flow that reacts (k mol/s), K t is 

the valve molar constant, x. is the molar fraction of species. 

Applying the Laplace transformation to the above equation 
and isolating the hydrogen partial pressure, we obtain the 
following expression: 

1 



— ^— (iV'V -2K r l) 



1 + r 



(2) 



where, T H2 - „ __, ,where, r„, expressed in seconds, 

is the time constant associated with hydrogen flow and is 
the function temperature, K is the constant dependent on 

Faraday's constant and number of electrons [Nj in the 
reaction. 



K, = 



N 
AF 



where p is the Faraday's constant. 



* **P»Fn*r- <-«»« 



^> 



.TlJ 



Fig. 2 SOFC Dynamic Model 

B. Thermal modeling 

The Nernst potential and loss mechanism are temperature 
sensitive, the temperature throughout the cell must be known 
to accurately determine the cell's electrical performance. 

In addition to the various types of heat transfer and the 
sources, the thermal model also includes the variation of 
material properties with temperature. The fuel cell power output 
is closely related to the temperature of the unit cell. The heat 
storage in the thin fuel unit or oxidant gas layer is neglected. 
The thin fuel unit or oxidant gas layers are lumped to the cell 
unit and gas layers are assumed to have the same temperature 
as the unit cell. 
E. Achenbach in [7] presents a relationship between two 

parameters S Q and F Q as given below 



F =0.725 



i.i 



(3) 



where ° 



ut Uj( l -v) ,„ 

and S = h. 



pC p h eg """ " tjAAT 
By simplifying the above eq. (3), we get 

!-. 0.909 mo c -1 

F Q = 0.72S Q 



■ff 



, 0.909 



= - 72 (%at)' 



(4) 



Where, 



uAi-ij) 

Y = — '- -h, 

7] 



if 



Uj is the power density (electric output), rj is the efficiency 

of the fuel cell, X is the thermal conductivity, Ar is change 

in temperature 

From the above eq. (4) we get, 

, 0.909 v 

(5) 



Ar-^ 



0.722 

This is the increase in temperature from the initial condition 
after a relaxation time t .Assuming the relation between 

t, dt and temperature j . We can predict the temperature at 
on next simulation instant as 



T Q =T 
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C. Nernst voltage equation 

Applying Nernst's equation and ohmic law (taking into 
account ohmic, concentration, and activation losses), the 
stack output voltage is represented by 



V J =V n -rI 



dc 



N, 



IF 



-In 



x H, x O, 2 



(7) 



(8) 



Where, (Vq) is open-circuit reversible potential (in volts), 
yE ) is standard reversible cell potential, (j^ J is mole 
fraction of species, (rj is ohmic resistances(in ohms), (i 7 ) is 
Faraday's constant (coulomb per kilo mole), [T ) is stack 
temperature in Kelvin, \N Q ) is the number of cells in 
stack, 77 , is activation losses in volts, 77 is concentration 

' I act ' I con 

losses in volts, yl) is the stack current in amperes. 

D. Power conditioning Unit Model 

The power conditioning unit is used to convert DC output 
voltage to AC. The power conditioning unit includes a DC/ 
DC converter to raise DC output voltage, followed by a DC/ 
AC inverter to convert the DC bus voltage to AC. In this 
paper, only a simple model of a DC/AC inverter is considered 
for the following reasons: the dynamic time constant of 
inverters is of the order of microseconds or milliseconds. 
The time constants for the reformer and stack are of the order 
of seconds. 

A simple model of the inverter is given in [9], where output 
voltage and output power are controlled using the inverter 
modulation index and the phase angle 8 of the AC 

voltage, V ac . The output voltage and the output power as a 
function of the modulation index and the phase angle can be 
written as: 

■mV ceU Zd 



V n 



P ac = ~ sin(<5J 



(9) 
(10) 



Q _ ( mV CCII f ~ m VccllV S COS (8) 

where, V ac is the AC output voltage of the inverter [ V] , m is 
the inverter modulation index, § is the phase angle of the 
AC voltage mV ce u [rad], P ac is the AC output power from 
the inverter [W], Q is the reactive output power from the 

inverter [VAR], V s is the load terminal voltage [V], x is the 
reactance of the line connecting the fuel cell to the infinite 
bus[Q]. 
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E. Active and Reactive Power Control 

In synchronous machines, the power angle is not 
measured, but the adjustment of the power angle occurs 
following changes in steam input and rotor speed. In the 
FCPP, there is no speed control but the similar relationship 
between output voltage phase angle and the flow of hydrogen 

can be adopted as follows. Given that the load voltage (VsZO) 
is constant and the AC source voltage out of the inverter V ac 
is given in (9), the angle § controls the power flow from the 
fuel cell to the load, as in (10). The phase angle can be 
controlled using the input flow of hydrogen. The expression 
for , therefore, provides the relationship between the power 
output as a regulated quantity, and the amount of flow of fuel 
input. This relationship is described by the following 
equations: 

Assuming a lossless inverter, we get 

Pac=Pdc=V cell I d2) 

According to the electrochemical relationships, a relation- 
ship between the stack current and the molar flow of hydro- 
gen can be written as 



Qh, 



N I 



2FU 
From equations (10), (12), and (13) 

. , 2FUX 

smo = q H 

mV s N 2 



(13) 



(14) 



Assuming a small phase angle sin 8 = 8 ', (14) can be written 

as 



2FUX 



<?#, 



(15) 
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Equation (15) describes the relationship between output 

voltage phage angle 8 and hydrogen flow </"#, . Equations 

(10) and (15) indicate that the active power as a function of 
the voltage phase angle can be controlled by controlling the 
amount of hydrogen flow. 

III. TS Fuzzy Controller 

The proportional plus integral (PI) controller requires 
precise linear mathematical models, which are difficult to 
obtain and may not give satisfactory performance under the 
transient conditions. The advantage of fuzzy logic controller 
is that it does not require a mathematical model of the system. 
Here in the paper, TS fuzzy controller is incorporated because 
of its simple structure. The linguistic rule consequent is made 

variable by means of its parameters. As the rule consequent 

is variable, the TS fuzzy control scheme can produce an infinite 
number of gain variation characteristics. In essence, the TS 
fuzzy controller is capable of offering more and better 
solutions to a wide variety of non-linear control problems. 
The deviations in the power error (error between reference 
and actual powers) are fuzzified using two input fuzzy sets P 
(positive) and N (negative). The membership function used 
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for the positive set is 



H P ( x , ) 





x i + L 
2L 



x i < - L 
- L < x, ■ < L 
x , > L 



Where, Xj (k) denotes the input to the fuzzy controller at the 
k th sampling instant given by 

x^k) = e(k) = P amf - P ac (17) 

and 

x 2 (k) = ^e(k) (18) 

For the negative set is 



f N ( * i ) = 



l x , < - L 

- x t + L 

2L 
x, > L 



L < x, < L 



(19) 




Fig. (3) Membership function for (a) Xy and (b) X 2 

The membership functions for and are shown in figures 
(3) (a) and (b) respectively. The values of L { and L, are chosen 
on the basis of maximum value of error and its integration. 
The TS fuzzy controller uses the four simplified rules as 

Rl: Ifxi(k) is P and x 2 {k) is P then 

ii\(k) = a\.x\{k) + a 2 .x 2 (k) 

R2: If x x Q€) is P and x 2 (k) is A/ then u 2 (k) = K 2 u x {k) 

R3: If xy(k) is^/and x 2 (k) is P then u 3 (k) = K 3 M,(fc) 

R4: If Xi(k) is N and x 2 (k) is N then u 4 (k) = K 4 u l (k) 
In the above rules, represents the consequent of the TS 
fuzzy controller. Using Zadeh's rule for AND operation and 
the general defuzzifier, the output of the TS fuzzy controller 

is. 



Z ( M j ) r « ;• ( * ) 



»(*) 



J-i 



Z ( " j ) 



where a = a^K, b = a 2 K and 

K = (Mi + K 2 /u 2 + K 3 /u 3 + ^ 4 ^ 4 )/(rt + /u 2 +/u 3 + /u 4 ) 

(16) ^ s ^' s operating condition dependent, the effective value 
of control gain varies widely during the control process. The 
values of a ; a 2 \ K 2 \ K 3 ; K 4 are given in APPENDIX-B. In the 
similar fashion TS fuzzy is incorporated to control reactive 
power. This concept is given in [10]. 

IV. Results and analysis 

Fuel cell is designed for 320 V D.C. output voltage. Fig. 4 
shows the output voltage of fuel cell. Temperature of fuel cell 
is dependent on active power requirement. Consider that 
active power is increased from 0.3 p.u. to 0.6 p.u at 1000 sec. 
and then decreased to 0.3 p.u at 2000sec. keeping reactive 
power constant; the temperature of fuel cell is changed 
accordingly as shown in Fig. 5. 

Fig. 6 (a) and 7 (b) shows the active and reactive powers 
which are dependent on each other. The active power 
requirement is changed from 0.2 p.u to 0.6 p.u at 200 sec. and 
again it is change d to 0.4 p.u at 300 sec. From Fig. 6, it is seen 
that fuel cell is capable of providing the active power 
requirement. As active power demand is changed, the reactive 
power will also change Fig. 7 shows the capability of fuel cell 
to meet the reactive power demand. In order to make reactive 
power requirement independent of active power, a decoupled 
control theory is used. Consider that the active power is kept 
constant at 0.5 p.u. as shown in Fig. 8 and reactive power 
demand is increased from 0.2 p.u to 0.6 p.u at 200 sec and 
then decreased to 0.4 p.u at 300 sec. Fig. 9 shows the capability 
of fuel cell to meet reactive power requirement which is 
controlled independently by applying decoupled control. 

Let initially rms value of active power is 30 KW and if the 
demand is that this rms value changes sinusoidally i.e., Pacref 
changes at 200 sees. Fig. 10 shows that the fuel cell with PI 
controller is tracking the required reference power (Pacref), 
but it does not track the reference power perfectly. 

To overcome this problem, TS fuzzy controller is 
incorporated (discussed earlier). Fig. 1 1 shows the tracking 
performance of fuel cell with TS fuzzy controller. From Fig. 
1 1 , it is seen that TS fuzzy controller is tracking the reference 
power perfectly. Now consider that reference active power 
changes (Pacref) shown in Fig. 12 and 13 shows the 
performance of fuel cell with PI controller and Fig. 13 shows 
the performance with TS fuzzy controller. From these figures, 
it is seen that TS fuzzy controller can meet the power demand 
perfectly as compare to PI controller. In the similar fashion 
the reactive power demand is changed and performance of 
fuel cell with PI controller and TS fuzzy controller is observed 
in figures 8.10 to 8. 12. From these figures, it is seen that TS 
fuzzy controller can track the required reactive power demand 
perfectly as compare to PI controller. 



(20) 



However for y =1, we get the centroid defuzzifier with u(k) 
given by 

u(k) = ax^(k) + bx 2 (k) 
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Fig. 4 Fuel cell output DC Voltage 




Fig. 5 Temperature of fuel cell 




Fig. 6 Active power 




Fig. 7 Reactive power 




Fig. 8 Active power (decoupled control) 




Fig. 9 Reactive power (decoupled control) 
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Fig. 10 Active Power with PI controller 
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Fig. 1 1 Active powers with Fuzzy controller 




Fig. 1 2 Reactive powers with PI controller 









Fig. 13 Reactive power with fuzzy controller 

V. Conclusions 

The proposed SOFC dynamic model is developed along 
with the subsystems of concentration block, activation block, 
water, oxygen and hydrogen block. Thermal block is also 
incorporated in the SOFC dynamic model. The integrated 
model includes fuel cell and power conditioning unit. The 
mathematical equations facilitate modeling DC to AC 
converter. The proposed model is tested with its active and 
reactive power outputs being compared with the required 
load demand. The dynamic behavior of SOFC model is 
analyzed by using PI and TS fuzzy controller separately. Here 
TS fuzzy logic controller is chosen because of its simple 
structure. It is observed that the dynamic performance of 
fuel cell with TS fuzzy controller is better in terms of tracking 
power demand as compared to the PI controller. 
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APPENDIX-A: Fuel Cell Operating Data 



Temperature[T) 


900 Z 'K 


No ioad voltage [£§ ) 


IV 


Number of cells [Nq) 


384 


Utilization factor [U ) 


0.8 


Line reactance ( X ) 


0.05 Q 


PI gain constants \K „ Kj J 


0.1 


Base Voltage {V r J 


400v 


Methane reference signal 


0.000015 kmois 


Base MVA 


100KVA 



APPENDIX-B 

For Active power control: o ; =10, a 2 ; =100, K 2 ;= K; =K=\ 
For Reactive power control: a =0.1, a' =10, K'= K:, =K=Q 
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